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Surface spontaneous parametric down-conversion is predicted as a consequence of continuity re- 
quirements for electric- and magnetic-field amplitudes at a discontinuity of \^ nonlinearity. A 
generalization of the usual two-photon spectral amplitude is suggested to describe this effect. Ex- 
amples of nonlinear layered structures and periodically-poled nonlinear crystals show that surface 
contributions to spontaneous down-conversion can be important. 



When studying the process of second-harmonic gen- 
eration under considerable phase mismatch more than 
thirty years ago, the generation of second-harmonic field 
from a boundary between two homogeneous media that 
differ by values of nonlinearity has been discovered 
0, 0] . The surface second-harmonic field arises here as 
a consequence of continuity requirements for projections 
of electric- and magnetic-field vector amplitudes into the 
plane of the boundary. Physically, a pumping field at 
frequency u> creates a step profile of nonlinear polariza- 
tion at frequency 2u> and with wave vector 2k{uj) that 
becomes the source of the usual volume second-harmonic 
field. The wave vector of the surface second-harmonic 
field is k(2u>) in agreement with dispersion properties of 
the nonlinear material. This effect is even found in non- 
linear media with negative index of refraction as the nu- 
merical solution of nonlinear Maxwell equations revealed 
in The studied parametric effect should be distin- 
guished from resonant surface second-harmonic genera- 
tion. 

Spontaneous parametric down-conversion (SPDC) m 
belongs together with second-harmonic generation to x' 2 ' 
processes. This poses the question about surface effects 
in SPDC. In volume SPDC, photon pairs are generated 
from the vacuum state, due to quantum fluctuations (or 
quantum noise) inherent in this state. In this case, a non- 
linear material responds to the presence of optical fields 
through quantum nonlinear polarization that acts as a 
source of new fields. In a close vicinity of the boundary, 
the interacting fields as well as the nonlinear polarization 
are modified in order to comply with natural fields' con- 
tinuity requirements at the boundary. This results in the 
generation of additional photon pairs from the area of 
the boundary (several wavelengths thick) that constitute 
surface SPDC. 

Our study of surface SPDC is organized as follows. 
Nonlinear Heisenberg equations are derived first to treat 
SPDC inside the nonlinear medium. Nonlinear correc- 
tions to electric- and magnetic-field amplitudes occur 
naturally at boundaries and give additional, i.e. surface, 
contributions to SPDC. Subsequently, the derivation of 
quantities characterizing the emitted photon pairs is ad- 
dressed. Finally, two important examples are discussed. 



Adopting the quantization of energy flux 0, Q we de- 
scribe the process of SPDC involving the signal, idler, 
and pump fields by the Heisenberg equations with an ap- 
propriate interaction momentum operator Gi n t [5[: 
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The positive- frequency part of an electric-field amplitude 
Em] can be expressed using annihilation operator a ma 
as follows (to = p,s,i; a = F, B): 
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0-rn a (z, UJ m ); 
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E m } = (iity. Subscript F (B) indicates a field prop- 
agating forward (backward), i.e. along +z (—z) axis. 
Symbol eo means permittivity of vacuum, d c g is effective 
nonlinear coefficient, A transverse area of the fields, c 
speed of light in vacuum, and h.c. replaces the hermitian- 
conjugated terms. Symbol k ma is a wave vector, ui m fre- 
quency, and n m index of refraction of field m a . 

The Heisenberg equations, e.g., for the signal- field op- 
erators a Sa (z,uj s ) can then be derived assuming equal- 
space commutation relations Q: 



da Sa (z,LO s ) 



dz 



E 



= ik Sa (u s )a Sa (z,ui s ) 

dwig{uj s , u!i)E^\0,uj s +u>i) 



0,-y=F,B' 

X exp[ik P0 ((J s + Ui)z]a\ (z,u)i), a = F,B; (3) 



Coupling constant g, g(u s ,uji) = 2id e sWu>s^>i / (cv 2tt 
\J n s (uj s )ni(uji)), is linearly proportional to nonlinear co- 
efficient d e g. 

The solution of Eq. |3|) for annihilation operators 
d Sa [z, u) s ) up to the first power of g gives us the formula 
for operator defined in Eq. @: 
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a Sa (0, u s ) + 



E 



du>ig(uJ s ,Ui) 



X-E^ (0, w s + Wi) &cp[iAkp jCn (u a , Ui)z/2] 
xzsmc[Akp^(ui s ,u} i )z/2}al_ i (0,u}i) ; a = F,B;(4) 

sinc(a;) = sin(x)/x and Afc ( g, Q7 (o; s , Wj) = k PlJ {uj s + u>i) — 
k Sa (u> s ) - ki y (ui). 

The positive-frequency magnetic-field amplitude op- 



erators H; 



can be derived using the formula 

—i/(tA} a fjio)9Ea^ (z,w a )/dz (fio denotes 
permeability of vacuum) provided that the electric-field 
[magnetic-field] amplitude E Sa [H Sa ] is polarized along 

+x [+y] axis. The obtained operator Hs„ can be de- 
composed into two parts denoted as H^^ and Hi„ nFl : 
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xE^\u> s + Wi) exp[ik Pj3 (uj a + Ui)z] 

x exp[-ifci 7 (wi)z]a^(0, a = F,B. (7) 



By definition, the magnetic-field amplitude operator 
,uj s ) is linearly proportional to the electric- 
field amplitude operator Ei^ (z, co s ). The remaining 

magnetic-field operator i?i^ nFr is of purely nonlinear ori- 
gin and the usual derivation of Fresnel relations does not 
take it into account. Standard approaches to nonlinear 
interactions thus do not involve this nonlinear term and 
so they neglect surface effects. We note that the 'nonlin- 
ear' magnetic-field operator H^ nFl occurs as a classical 
field amplitude also in the description of stimulated para- 
metric processes (e.g., in difference- frequency generation) 
and yields surface contributions to these processes. 

The electric- and magnetic-field amplitudes 
E ma (z, uj rn ) and H m<x (z, uj m ) originating in the nonlinear 
interaction and written in Eqs. 0$ and ([SJ— [7]) have to 
obey continuity requirements at the input and output 
boundaries of the nonlinear medium. We illustrate our 
approach to this problem considering the signal field at 
the input boundary (z = 0). Four electric and magnetic 
fields are involved in the continuity requirements at this 
boundary (see Fig. [1]): two at the linear left-hand side 
[denoted by superscript (0)] and two at the nonlinear 
right-hand side. Because the magnetic-field amplitudes 
H~ and H, inside the nonlinear medium have also 
nonlinear contributions -ff"J r and -ff"J r given in Eq. Q 
additional (surface) amplitude corrections SE SF and 
8E s °b [together with SH SF and SH aB '] in the fields 
leaving the boundary naturally occur. The amplitude 
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FIG. 1: Scheme showing electric and magnetic fields at the 
input boundary. For details, see the text. 



corrections SE^ and 5Hs B of the outgoing field outside 
the nonlinear medium can be involved in the fields obey- 
ing Fresnel relations [1| at the expense of introduction 
of fictitious amplitude corrections SE SB and SH SB of 
the field impinging at the boundary from its nonlinear 
side. A detailed analysis then results in two equations 
for the surface amplitude corrections of fields inside the 
nonlinear medium: 

= 5E SF (p)-5E SB (0), 

= iO0) + SH tp (0) + H^(0) SH SB (0). (8) 

The positive-frequency parts of surface amplitude- 
correction operators SE^ and 5H a ^ occurring in the 
quantum form of Eqs. 1(5)) can be expressed using 
annihilation-operator corrections 5a Sa similarly to the 
corresponding amplitude operators E^ and H^ in 
Eqs. © and ©. The solution of Eqs. © for Sa SF and 
5a SB then takes the form: 



5a sp (0,u> s ) = <5a SB (0,w s ) = 



(kJi 9 {uj s , u>i)E$ (0, oj s + Ui)al (0, oj,). (9) 



Similar considerations appropriate for the output 
boundary leaves us finally with an expression for opera- 
tors a Sa (L,uj s ) valid up to the first power of g (L stands 
for the length of nonlinear medium): 



0,1=F,B 



d^i Fl,fa{us,Ui)a£** (L,tJi), a = F,B 



(10) 



Operators a, (L, u> s ) correspond to free-field linear 
propagation, i.e. without photon-pair generation. The 
idler-field amplitudes can be analyzed along the same 
vein. 

The generalized two-photon spectral amplitudes T s 
and J- 1 defined in Eq. (fTOf describe properties of a gener- 
ated photon pair and are composed of two contributions: 



(ii) 
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Two-photon spectral amplitude .F™ 1 of the volume con- 
tribution has the well-known form: 

x exp[ik Pos (uj s + uJi)L] exp[— iAk a ^ 7 (w s , u>i)L/2] 
xLsinc[Afc Qi/ 3 7 (t i ; s ,w l )L/2]; a,0,~/ = F,B. (12) 

On the other hand, surface contributions J^ m ' su " to the 
two-photon spectral amplitudes can be expressed as: 

= v^K,^)^,^,^), (is) 

where 

V™,i3y{UB,Ui) = 7^7 ; ; a,p,j = F,B. (14) 

"to J 

The structure of surface contributions as described by 
the two-photon amplitudes jF 5 > sm ' f and jF J > surf resembles 
that of the volume contribution as the formula in Eq. (|13p 
indicates. The physical interpretation is as follows. At a 
boundary, the only restriction for photon-pair generation 
is imposed by the conservation of energy. However, the 
mutual interference of two-photon amplitudes originating 
at the input and output boundaries leads to the result 
that resembles the usual phase-matching conditions. We 
note that lim L ^ ^* m ' sur£ = 0. 

We further consider photon pairs with both photons 
propagating forward and use operators a m {u> m ) (m = 
s,i) defined outside the nonlinear medium. The joint 
signal- idler photon- number density n(u> s ,L0i) at the out- 
put plane of the nonlinear medium is given as: 

n(iJ s ,Ui) = ^ a\ (oJ s )a s (a; s )aj (u^a* (wj) +h.c. ^ /2. 

(15) 

Symbol ( ) denotes averaging over the initial signal- and 
idler-field vacuum state. Introducing two-photon spec- 
tral amplitudes T s and T l (transmission coefficients t m 
describe the output boundary), 



(16) 



T m (uj s ,uJi) = t s (uj s )ti(uji)J 7 F l FF (uj s ,Lu l ), m = s,i, 

we arrive at the following formula: 

n(uj s ,cJi) = Re{J rS *(o; s ,LJi)J : " 4 ( w s,^i)}. (17) 

As this example illustrates, a generalization of the usual 
formalism based on a two-photon spectral amplitude can 
be given providing formulas for all physical quantities 
characterizing photon pairs. 

The volume interaction among the forward- 
propagating pump, signal, and idler fields dominates in 
bulk nonlinear crystals several mm long. According to 
our model, the surface contributions can be approxi- 
mately included into the usual formalism working with 
a two-photon spectral amplitude $ vo1 (see, e.g., [3, Il0| ) 
using the formal substitution: 
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FIG. 2: Signal-field spectra S™ 1 (solid curve denoted as a) 
and ST 1 (solid curve denoted as b) of volume and surface 



SPDC, respectively, and ratio S, 



vol+surf 



/ Ss° of the spectra 



xJl + V 



F,FF 



K,^)$ vol K,^:); (is) 



with (5T 1+surf ) and without (5J o1 ) the inclusion of surface 
SPDC (dashed curve); S s (us) = huj a J du)i n(io B ,LOi). 



V is defined in Eq. fH]) . If the nonlinear interaction is 
perfectly phase-matched {AkF,FF(u° , cj°) = 0] the sur- 
face contributions at central frequencies are zero. 

Contrary to the bulk nonlinear crystals, surface SPDC 
cannot be neglected in nonlinear layered structures com- 
posed of layers typically several hundreds of nm long. In 
this case all possible nonlinear interactions as described 
by the momentum operator G m t in Eq. (fT]) give appre- 
ciable contributions. Fulfilment of phase- matching condi- 
tions is not important here, because Akl -C n (I denotes 
a typical length of one layer). Surface SPDC similarly as 
volume SPDC from an individual layer is weak but both 
of them arc highly enhanced by constructive interference 
of fields from different layers. A generalization of the 
presented theory to layered structures is straightforward 
following the work presented in [ll|, [l2j • 

As an example, we consider a structure composed of 25 
layers of nonlinear GaN 117 nm thick that sandwich 24 
linear layers of A1N 180 nm thick and studied previously 
in [ll|. Volume SPDC gives efficient photon-pair gener- 
ation at degenerate signal- and idler-field frequencies for 
the signal- field emission angle 14 deg [ll[ (see Fig. [2]) as- 
suming a normally incident pump field at X p = 664.5 nm 
and s-polarized fields. Additional photon pairs originate 
in surface SPDC. Their intensity is cca 20 % of that com- 
ing from the volume. However, both contributions are 
in phase and add constructively so that the inclusion of 
surface SPDC roughly doubles the number of emitted 
photon pairs in the spectral area of efficient photon-pair 
generation (see Fig. [2]). Different contributions to surface 
SPDC can be quantified using coefficients V^ 7 defined 
in Eq. (TT4")) . Whenever the lengths of nonlinear layers are 
less or comparable to the coherence length of the non- 
linear process, we observe appreciable contributions of 
surface terms. For example, the coherence length equals 
approximately 1 /im for GaN in our case. 

Surface effects give also an important contribution to 
photon-pair generation rates in periodically-poled non- 
linear materials with sufficiently short poling periods. 
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FIG. 3: Relative contribution n vo1+sut1 /N vo1 - 1 of surface 
SPDC to the photon-pair generation rate (points •) and in- 
verse A'j 1 of the optimum poling period giving quasi-phase- 
matching for A° = A- = 2X° (solid curve) as they depend 
on cw pump-field wavelength Ap 1 in LiNbOs 5 mm long; 
N = J dui s J du>i n(u s ,Ui). 

Here, as an example, we consider frequency-degenerate 
SPDC in periodically-poled L1NDO3 with the optical axis 
perpendicular to the direction of collinearly-propagating 
fields; their polarizations are parallel to the optical axis. 
Whereas the surface effects contribute to photon-pair 
generation rate _/v vol + surf only by several percent for the 
pump wavelength Xz, = 1 /zm, the increase of photon-pair 
generation rate TV by 50 % is observed for A° = 0.35 /jm 
(see Fig. As the curves in Fig. [3] indicate the rela- 
tive contribution 7V vol+surf /A vo1 — 1 of surface terms is 
roughly proportional to the inverse 1 / A n i of poling period 
that is linearly proportional to the density of surfaces per 
a unit length. We note that domains shorter than 1 /im 



can be fabricated using light-induced domain engineering 

m, _ 

Surface SPDC occurs also in nonlinear wave-guiding 
structures, i.e. under the condition of total reflection and 
presence of evanescent waves. A detailed analysis has 
shown that the formulas presented above remain valid 
also in this case provided that we consider propagation 
constants (3 instead of wave vectors. Coupling constant 
g then involves the overlap integral over mode functions 
of the interacting fields. This is particularly interesting 
for nonlinear photonic-band-gap fibers. 

Surface SPDC is by no means restricted to ID non- 
linear structures: even greater relative contributions are 
expected in 2D and 3D nonlinear samples. Surface ef- 
fects will also affect stimulated x processes like second- 
harmonic or second-subharmonic generation when stud- 
ied under comparable conditions. Qualitatively, they will 
effectively enhance the nonlinearity. This may be partic- 
ularly interesting for squeezed-light generation. 

In conclusion, surface SPDC has been predicted. Gen- 
eralized signal- and idler-field two-photon spectral am- 
plitudes have been suggested to determine properties 
of emitted photon pairs. Surface SPDC is important 
whenever strongly phase-mismatched nonlinear interac- 
tions give considerable contributions. This occurs, e.g., 
in nonlinear layered structures or periodically-poled ma- 
terials where surface and volume contributions can be 
comparable. Surface SPDC may affect optimum design 
of these structures that are considered as promising ver- 
satile sources of photon pairs for optoelectronics. 
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